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ABSTRACT 

We present a new analytic estimate for the energy required to create a constant density core within 
a dark matter halo. Our new estimate, based on more realistic assumptions, leads to a required 
energy that is orders of magnitude lower than is claimed in earlier work. We define a core size based 
on the logarithmic slope of the dark matter density profile so that it is insensitive to the functional 
form used to fit observed data. The energy required to form a core depends sensitively on the radial 
scale over which dark matter within the cusp is redistributed within the halo. Simulations indicate 
that within a region of comparable size to the active star forming regions of the central galaxy that 
inhabits the halo, dark matter particles have their orbits radially increased by a factor of 2-3 during 
core formation. Thus the inner properties of the dark matter halo, such as halo concentration, and 
final core size, set the energy requirements. As a result, the energy cost increases slowly with halo 
mass as for core sizes < Ikpc. We use the expected star formation history for a given dark 

matter halo mass to predict dwarf galaxy core sizes. We find that supernovae alone would create well 
over 4kpc cores in 10^° Mq dwarf galaxies if 100% of the energy were transferred to dark matter 
particle orbits. We can directly constrain the efficiency factor by studying galaxies with known stellar 
content and core size, such as Fornax. We find that the efficiency of coupling between stellar feedback 
and dark matter orbital energy need only be at the 1% level or less to explain Fornax’s 1 kpc core. 
Subject headings: dark matter — galaxies: dwarf — galaxies: evolution — galaxies: formation — 
galaxies: halos — Local Group 


1. INTRODUCTION 

Collisionless simulations of dark matter halos in 
ACDM cosmology consistently predict a den- 


Dubinski & Carlberd 

19911: INavarro et al.l 119951 119971: 

Bullock et al.1 120011: 

Klvoin et al.l 120011: IStadel et al.1 


rolls over to a steep slope at large radii where it is 
broadly consistent wit h the rotation curves of man y 
late-type galaxies (e.g. iRubin et al.lIlQSDl: iBosmalllQSlD. 
Howe ver, the inner cusp (e.g. p oc r“^, iNavarro et al.l 
is inconsistent with observations of rotation 


profiles in dwarf galaxies both within the local group 
and beyond (e .p. Flores.^i^Primac2jl99^_lMoor 
BurkertI 


Cote et al.l 120001 : 

[2007 


iGilmore et al 

lOh et al. 20111: 


Kuzio de Narav et al 


_ .Kuzio de Narav et ah. 

Walker &: Penarrubial 120111 : 


Amorisco et al.ll2013l : lAdams et al.ll2014li . Instead, these 
dwarf galaxies are better fit using a constant density 
within the central kpc or so. Two classes of solutions 
have been proposed to resolve this discrepancy. The 
first suggests that coupling to baryon s can reshape 
the da r k matter cusp into a core (e.g. INavarro et akl 
1996at [El-Zant et al.l l200il iWeinberg fc Katd ^0^ 


Read fc Cilmorel 1200,^ iMashchenko et all 1200611 . 

These solutions have met w i th some criticism (e.g . 
Sellwoodll2003l: lTasitsiomil[2003l : lJardel fc Sellwoodll2009l: 


Dubinski et al.l1T009ll . In particular. iPefiarrubia et ^ 
MTl " claimed it was not energetically feasible. How¬ 
ever, the sub-structure present within massive halos 
can induce collisionless scattering of dark matter 
particles within the centre, flattening the inner cusp 

aj max@mcmaster. ca 


witho ut the need of baryons (e.g. iMa fc Bovlan-KolchinI 
l2004|l . Alternative models to ACDM have been also 
been suggested, such as ‘warm’ dark matter (e.g. 
iHogan fc Dalcant^ I2000I1 or (no n-gravitational) dark 
matt er self-interactions (e.g. iSoergel fc Steinhardd 
[ 2 (^. 

IMashchenko et all (|2006ll argued that cyclic star 
formation bursts within the centres of dark matter halos 
could provide the necessary energy to re-shape the dark 
matter cusp into a dark matter core. In this model, 
feedback from type H supernovae (SNe) periodically 
drives bulk gas motions in the star forming regions in the 
central kpc of dwarf galaxies. Since gas can dominate 
the central gravitational field, these bulk motions can 
significantly perturb the central gravitational potential 
on timescales comparable to dark matter orbital times 
leading to an efficient transfer of kinetic energy to the 
central dark matter distribution and the formation of a 
significant dark m a tter c ore. A notable feature of the 
IMashchenko et all (|20n6[ i framework was that the gas 
need not be expelled or even pushed very far outside the 
core-forming reg ion (as opposed to the impulsive blow¬ 
out s c heme e.g. [Navarro et al. Il996al: iRead fc Gilmord 
120051) . iPontzen fc Covernatol ( ^121) gave an analytical 
description of how changing the dominate potential 
leads to irreversible heating of the dark matter orbits. 

The mechan ism was first verif ie d in a cosmological 
simulation by IMashchenko et all (|2008f) . A required 
feature for simulations to flatten cusps is realistically 
clustered star formation so as to ensure a star for¬ 
mation history that significantly fluctuates in both 
space and time. Several simulations since have repro¬ 
duced these results over a range of simulated galaxy 
masses (e.g. iGovernato et al.l[2010L 120121: iZolotov et al.l 
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20121: ITevssier et alJ 

2013; 

Pontzen & Governatol 

2014) 


The process also has 
consequences for other collisionless components such as 
stars and star clusters. Initially centrally concentrated 
populations (such as younger stars) are pushed outward 
over time, leading to older populations having larger 
orbits on average. This results in a broad match to 
observations of the s t ellar content of dwarf galaxies (e.g. 
iMaxwell et al.l l2012t iTevssier et al.l 120131 : ?). It also 
has interesting consequen ces for multiple popu lations in 
g alactic globular clu s ters (Maxwell et al.l 120141 1 . 

iPefiarrubia et al] (j2012 1 made the first analytical 
estimate of the energy difference between cored and 
cuspy dark matter halos of the same mass, finding that 
SNe alone would have difficulty providing the energy 
required to form dark matter cores given the average 
stellar masses of dwarf galaxies. Though there are other 
sources of energy available from star formation, such as 
stellar winds and the ionizing flux from young stars (e.g. 
iLeitherer et al.l 1199911 . it is not clear that this added 
energy would be sufficient to remove the discrepancy or 
couple as well as SNe-driven bulk gas motions. How¬ 
ever, the IPefiarrubia et al.l (j2012l l results are difficult 
to reconcile with the core forming simulations which 
typically employ feedback only from SNe (and with 
realistic efficiencies). 

In this work we re-examine the energy requirements 
and show that the formation of dark matter cores is en¬ 
ergetically plausible in a ACDM cosmology. The layout 
of the paper is as follows. |J2] deals with the ambiguity 
inherent in the definition of a core size, and argues for 
a new, general way to define the core radius. In SJ3] we 
discuss the energy required to form dark matter cores in 
halos that were initially cuspy using assumptions that 
are more consistent with our understanding of how the 
process operates. In ij4] we show predictions for core 
sizes in dark matter halos and examine the implications 
for observed systems such as Fornax. 


2. CONSISTENT DEFINITION OF CORE SIZE 

Many observed cores in dwarf galaxies are in¬ 
ferred from the rotation profiles of the stars and 
gas by fin ding the best-fit dark - matter density 


model (e.g. iFlores fc PrimykI 11994 iBurkerd 
Kuzio_deNara^_e£^ 20061: IG ilmore et al.l 


1995 


2007 


Kuzio de Narav et al.l 120081 : lOhet al.ll2011[l . Commonly. 


one would use either a cored isothermal sphere, whose 
density prof ile goes as at large radii, or that 
proposed by iBurkertI (1199511 where the density goes as 
r~^ at large radii, 


p{r) oc 


1 


(1 + (*))(! +(*>4 ■ 


( 1 ) 


Both satisfy the requirement that the density profile 
transitions to a flat core within a scale radius. Without 
probing the rotation profile to large radii it is difficult 
to determine which density profile is the best fit. The 
cored isothermal pr ofile predicts a constant velocity at 
large radii while the IBurkertI (119951 1 profile looks like th e 
standard galaxy rotation curve (e.g. iRubin et al.lll978ll . 
However, the common practice of defining the core size 
as the scale radius leads to large variations in core radii 
dependent on the chosen functional form. 


To avoid this ambiguity, we define the core radius as 
the radius, r^, at which the logarithmic slope of the den¬ 
sity profile satisfies: 

d Inp 1 

dlnr 2 


We will use this definition from now on when referring 
to the core size. Characterizing cores this way is con¬ 
sistent with the density profil e slopes found in simu¬ 
lations of dwarf galax i es ('e.g. iMa shchen ko et al.l 120081 : 
IGovernato et al.l 120101 : ITevssier et al l 1201311 and with 
the mass profile slopes observed in local group dwarfs 
le.g. IWalker fc Peharrubi^ 120111 : lAmorisco et al.l 120131 : 

I Adams et al.ll2014ll . Furthermore, the simulated galaxies 
always transition to a density profile with a strongly neg¬ 
ative slope (din p/d In r ~ —3) beyond a few kpc where 
the effects of baryonic feedback are limited. This defi¬ 
nition of the core radius lessens the reliance on an as¬ 
sumed dark matter density distribution, and it can be 
determined directly from the velocity profile. It is also 
straightforward to relate r^ (for any given density profile) 
to re- 

To demonstrate the robustness of our new core defini¬ 
tion, and to illustrate the ambiguity with using the scale 
radius, we show in Figure [T] the rotational velocity pro¬ 
file for th e Low Surface Brightness (L SB) galaxy NGC 
959 from iKuzio de Narav et al.l (1200811 . These authors 
assumed a cored isothermal density profile, and the best 
fit of this profile to the data is shown as the solid line. 
Although it is perfectly sufficient to match the profile 
in the inner 1-2 kpc, it does not roll over sufficiently to 
match the slope in the outer part of a cosmological halo. 
We have also fit three other density profiles. The first, 
show n as the dotted line in the figure, is (e.g. iWidrowl 
l2?ionL ’Fixed’ in the figure): 


p(r) oc 


1 

( 1 + (*)")“'"■ 


( 2 ) 


The second, shown as the dashed line, is a pseudo- 
isothermal profile, where we let the power index vary: 


pW « 1 , Lu - (3) 

v rs / 

Lastly, we show the lBurkeiil (|1995l Eqn. [T]) profile as the 
dash-dot line. 

The differences between the three rotation profile fits 
are only evident at radii past the furthest observable 
point, where the predicted rotation velocity rolls over. 
We have also done th i s for the eight other LSBs from 
IKuzio de Narav et al.l (1200811 (see the Appendix), and 
the resulting scale radii and core radii are shown in Table 
[TJ It is evident that velocity rotation data alone are not 
sufficient to distinguish the best fit dark matter density 
profile. Each fitting function requires a different scale 
radius, the common definition of the core size. However, 
using d Inp/d Inr to define the core radius leads to consis¬ 
tent core sizes, independent of which of the four density 
profiles was used, so long as the core location is comfort¬ 
ably within the range of data points. _ 

The re are some LSBs in the IKuzio de Narav et al.l 
(|2008ll sample that do not show consistent Tc values, par¬ 
ticularly UGC 4325 and DDO 64. These are galaxies 
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where the core size is beyond the edge of the data, so a 
fitting algorithm is only able to place lower bounds on 
Tc and is certainly not able to constrain the parameters 
of a non-linear density function. 


3. THE CORE FORMATION ENERGY BUDGET 

In this section we estimate the energy required to 
build a core in an initially cuspy density profile. We 
assu me that each cuspy halo begins w ith an ‘NEW’ pro¬ 
file (jNavarro et al.lllQ^ ll996bL Il997fl . If we assume the 
halos are virialized, we can estimate the energy as the 
difference between the work required to form the cusp 
and the work required to form the core, as was done in 
iPeharrubia et al.f (|2012ll : 


AE = 


2 


(4) 


where Wh and Wc are the cusp and core potential en¬ 
ergies, respectively. AE > 0 implies the cusped dark 
matter halo had to gain energy to redistribute mass to 
form a core. If the simplifying assumption is made that 
both halos are spherically s ymmetric, we can write th e 
potential energy integral as (|Binnev &: Tremainell200^ : 


and reflect the collisionless formation process. This pic¬ 


ture i s supported by simulations (e.g. iMashchenko et al 


2006 . 2008HGoverna.to et a.l.ll2010l:lPontzen fc Governato 


2012 : Tevssier et al.l 201^ Mada u et al.l 2014 ). Previous 


simulations ( Mashchenko et al.ll2008HMadau et al.ll201^ 

suggest that for a ^ 10^° Mq halo Tm is about 2-3 times 
the core radius. 

Our method proceeds as follows: 


• Since we do not have any robust predictions of what 
the core profile should be, we adopt Equation © as 
the density profile when computing Wc- (We have 
verified that using a different density profile, such 
as Equation o or Equation m, does not change 
our results.) 

• We use the relations g i ven b y IMaccio et al.l (120071) 
and IBrvan &: NormanI (1199811 to set the NEW pro¬ 
file parameters at redshi ft zero as a function o f 
virial halo mass, Mh (see IPeharrubia et al.ll20O) . 
We then have p{rni) and M{r^) as functions of M^. 

• We find the set of parameters po, S, and rc for 
the cored density profile that simultaneously satisfy 
pi^m) and M(r-[a) from the previous step. 


W = -47rG 


rp{r)M(r) dr. 


0 


(5) 


For a density profile that goes as p oc r~^ at large r, the 
integral converges for /3 > 5/2; all of the cored profiles 
considered (if 5 > 5/2 in Eqn. [3]) and the NEW profile 
lead to well defined values for W. 

Taking the upper limit in Equation ([5]) to infinity leads 
to a very large estimate for the energy (quite apart from 
the ambiguity of embedding the halo in the cosmic back¬ 
ground) and is not appropriate given that we are con¬ 
sidering redistrib uting mass withi n the central part of 
an existing halo. IPeharrubia et al] ()2ni2[ l chose the halo 
virial radius rh as the upper limit for the potential energy 
integral, and normalised the mass of their cored density 
profile, 


p(r) 


{rc + r){rs + r)^ ’ 


( 6 ) 


to that of their initial NEW halo at rh- However, simply 
normalising both halos to have the same mass at rh does 
not remove the requirement that Equation ([S]) be carried 
to infinity. In order for the potential energy integral to 
be valid with a finite upper limit r^, the density and 
mass profiles for both halos must ma tch at that radius. 
The choice of iPeharrubia et all (|2012f) to match the mass 
at rh and to truncate the integral for W at that radius 
builds in an implicit discontinuity in the density at that 
radius. Furthermore, moving the cusp mass to near the 
virial radius of the halo leads to a core formation energy 
budget from Equations ([S]) and dU which will be signifi¬ 
cantly overestimated compared with the more restricted 
redistribution considered here (r„i<C rh). 

The motivation for constraining the mass and density 
profiles of the pre- and post-core halos outside rm to be 
the same is simply that we expect that strong coupling 
of star formation feedback to dark matter is only effec¬ 
tive in, and only modifies, the central parts of the halo. 
The outer parts of the halo will remain largely untouched 


• We set rju = 3rc, where r^ is given by our new 
definition. Figure [3] illustrates the profile and mass 
matching. 

• We compute the potential energy difference in 
Eqn dH). 

The preceding steps give AE in terms of rc and Mh. 


4. RESULTS 

In Figure [3] we show our calculation for the amount 
of energy required to turn a dark matter cusp into a 
core as a function of halo mass. Since we may write 
approximately: 


Wc 


GM{< r^f 


(7) 


for the cored profile, we see, as expected, that the pa¬ 
rameter that most influences the energy required to con¬ 
vert a cusp into a core is rm oc rc. A larger core moves 
more mass to larger radii and thus requires more energy 
to redistribute it from the cusp. On the other hand, 
if we allow rm to scale as a fixed fraction of rh (we 
chose rm/rh=l/6), we find that AE scales as Mh^ ®. 
This scaling is shown by the solid black line in Figure 
51 and lies within the r ange of solutions calculated by 
Peharrubia et al.l d2012l l (the gray shaded region in the 
figure). A more realistic estimate for rm, as a small mul¬ 
tiple of rc (here 3), leads to a greatly reduced energy 
requirement and a much shallower dependence on halo 
mass. This is shown by the grey lines in Figure [3] for 
several core sizes. Our calculations suggest the energy 
scaling with halo mass is much weaker when rc (and 
hence rm) is fixed. In the 10®-10^° Mq halo mass range, 
AE scales as Mh°'^ for a 1 kpc core and flattens to Mh°'^ 
for a 100 pc core. This flattening is related to how the 
mass interior to rm depends on the halo concentration. 
In the functional form we adopted the concentration is 
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r [kpc] r [kpc] 

Fig. 1.— LEFT: The velocity profile of the LSB NGC 959 from IKuzio de Narav et al.l H2008li . with the best fit ‘ zero disk’ cored isothermal 
profile shown as the solid line. The dotted line shows the best fit using Equation m, the dashed line shows Equation and the dash-dot 
line shows the best fit using the IBurkerU II1995I) profile (Eqn. m- For many of these halos, differences between the cored profiles are only 
evident at radii well outside the available observations. RIGHT: The dark matter density profiles derived from the best fit velocity profiles. 
For the same observed dwarf galaxy four different scale radii, the traditional definition of the core size, can be derived. The filled symbols 
show Vs for each of the four profiles, while the open symbols show Vc, where dlnp/dlnr = —1/2. 


TABLE 1 

Best fit parameters. 




Cored 


Pseudo 



Fixed 


Burkert 

Galaxy 

S 

Vs 

Tc 

5 

Ts 

Tc 

a 

rs 

rc 

Ts 

rc 

NGC 959 

2 

0.4 

0.23 

^•-LG_o.i8 

n 42+0.02 

G.9:Z_q Q2 

0.24 

i.6ol°:l? 

0.68t°:°5 

0.25 

0 74+0-01 

G. i O-o.oi 

0.27 

NGC 7137 

2 

0.6 

0.35 


0 61 +0 08 
G.Ol_ 0 .i 0 

0.35 

l-4llo.|6 

l-14l°1^ 

0.36 

1 qq+O.04 

-L.uy_Q Q4 

0.40 

UGC 11820 

2 

1.1 

0.64 


1 1 c:+0.06 

-l--iO_o 06 

0.70 

2.76^" 

0.89t°:gi 

0.50 

9-11 +0.28 

^•-*--*--0.23 

0.77 

UGC 128 

2 

2.3 

1.33 

1 Q7+0.36 
-L.y < -0.32 

2.48t°;|? 

1.43 

1 Q9+0-37 
-L.oz,_o 27 

5 01 +108 
o.Ui_o.94 

1.48 

4.50+;l® 

1.64 

UGC 191 

2 

1.7 

0.98 

-L-^^_0.64 

1.84t°;«5 

1.11 

-L-^o_0.48 

q q'T+6.26 

0 . 0 ^ _i 67 

1.14 

q rq+ 0.68 
'^•'->o_o 50 

1.28 

UGC 1551 

2 

1.3 

0.75 

-L.oy_Q 33 

1 qo+O.lS 

-L.o^_0.3i 

0.77 

1 4q+i-45 

-L-9:y_o 49 

2 oq+2.45 
z,.oy_o.9o 

0.81 

2.40+;?i 

0.87 

UGC 4325 

2 

4.1 

2.37 

1.74-0.35 

2.94==- 

1.74 

2.05-0.30 

2.52® 

1.08 

9.52® 

3.46 

DDO 64 

2 

3.3 

1.91 

9.51-7,89 

7.44-5.52 

5.49 

8.83-6.98 

7.44-5.77 

6.20 

8.76® 

3.18 

F583-1 

2 

2.5 

1.44 

2 1(7+0.69 
^•-L'J_0.53 

9 crr + O.ie 
^•O0_o.i9 

1.46 

1 07 + 0.92 
-0.55 

q cry+l.OS 
0.0 1 _0.77 

1.51 

4.49+;1«5 

1.63 


Note. — These parameters represent the best fitting density profi les required to match the central core density out 
to the edge of the rotation profile data of IKuzio de Narav et ] ll200l^ . All radii are in kpc. The scale sizes for the four 
different profiles, the traditional definition of the core size, vary widely despite having the same inner central density 
and matching the velocity data. The new definition of core size is much more consistent among profiles. Error bars 
are determined from least squares minimization of the model velocity profile residuals (see Appendix). Since the fitting 
algorithm uses rg and the power law index as parameters, it is difficult to put error bars on rc (which is derived directly 
from the logarthimic slope). However, the errors on Vc should be similar to those on r^. 

^ Lower limits. 


inversely proportional to Mh l|Maccio et al.ll200^ . Vary¬ 
ing the ratio r^jr^ creates a modest vertical shift in the 
Ai? curves. 

In Figure [3] we also plot the stellar ma ss to halo 
mass (Ms-Mh) relatio ns from iKravtsovI (1201011 fdashed) , 
iBehroozi et al.l (120101) fdotted!. and iMoster et al.l (|2010[! 
(dash-dot) at z = 0. W e have converted Mg to total SNe 
energy using a iKroiioal (l200l IMF, assuming that every 
star above 8 Mq contributes 10®^ ergs to the energy bud¬ 
get. Our results show that there is no confliet between 
the amount of energy available to form a dark matter 
core and the ACDM framework, regardless of halo mass. 
The main question is how large the dark matter core is 
likely to be. 

In reality, not all of the SNe energy will couple to the 
cusp dark matter through the induced perturbations to 
the gravitational potential. This inefficiency can be rep¬ 
resented by scaling the Mg-Mh relations shown in Fig¬ 
ure |3] by some coupling efficiency factor, e, (for example, 


the 40% value used bv iPefiarrubia et aI1l2012l! . Only for 
e < 0.01% would there be no possibility of creating a core 
larger than 250 pc in the 10^-10 ^'’Mq halo range. The 
implication is that most, if not all, dwarf galaxies should 
host dark matter cores if their central star formation rate 
was high enough to generate 10^-10^ Mq in stars. 


4.1. Efficiency of Core Formation 

In the previous section we showed that stellar feed¬ 
back, of which we have only considered SNe, can 
provide enough energy to the dense star forming gas 
in the centres of dwarf galaxies to transform the dark 
matter cusps predicted by numerical experiments to the 
flatte ned dark matter co r es inf er red from observations 
(e.g. iMashchenko et aP I2008I iPontzen &: Governatol 
l2012t iMaxwell et al.l 201211. This is contrary to previous 
work le.g. IPeharrubia et al.l l2012fl and results from 
our adoption of a more consistent mass normalisation 
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Fig. 2.— A pictorial representation of our procedure to match 
both density and mass of the core and cusp density profiles at rm- 
The leftmost shaded region shows the amount of mass that must 
be redistributed from the cusp to larger radii (up to rm) in the 
cored profile. Not having to redistribute mass to larger radii near 
the virial radius explains why our predictions for the energy cost 
of core formation are significantly below previous estimates. It is 
purely coincidental that in this illustration the core radius, r^, is 
close to the radius at which the two density profiles first cross. 


criterion. This was motivated in part by the results of 
hydrodynamic simulations and by the observation that 
the most efficient injection of star formation feedback 
energy will occur at the deepest point in the galactic 
potential, which implies a physical limit to the radius 
to which cusp dark matter can be redistributed. In our 
calculations, we have parameterized this limit as the 
radius rm. Our results suggest significant dark matter 
cores can be formed with less than 1% of the total 
SNe energy available over the star forming history of a 
typical dwarf galaxy. 

So far we have been concerned with whether a dark 
matter core can form in an initially cuspy halo due to 
star formation feedback. Our new calculations allow us 
to ask a more interesting question: what is the typical 
core size for a given halo mass? In Figure |4] we show the 
predicted core size as a function of halo mass using the 
three Mg-Mh relations shown in Figure [3] and 100% of 
the associated SNe energy. Although Figure |4] reinforces 
our main conclusion that dark matter cores can form 
from cusps in even the smallest dwarf galaxy halos, 
it also emphasizes the uncertainty in trying to make 
detailed predictions. The variation in t he low mass 
slope of the rel a tion b etween the iMos ter et all 

(|2()10D . iBehroozi et al.l (|2010t ). and lKravtsoX ~( 2010f ) 
results corresponds directly to a factor of 2-4 variation 
in predicted Tc as a function of halo mass. Furthermore, 
these relations only give the average stellar mass as a 
function of halo mass. 

Figure m is further complicated because we have not 
yet considered the efficiency with which the energy from 
stellar feedback is transferred to dark matter particle 
orbits. The coupling efficiency will likely vary in time 
with changing conditions in the centre of the halo, and 
will depend on the star formation rate, gas accretion 
rate, gas heating rate, and t he dark matter particle 
orbits; the value of e used in iPenarriibia et ahl (I2012D 
is merely a time-averaged value over the dwarf galaxy 
star formation and gas accretion histories. Finally, our 
results show that the energy required to form the core, 
AE, has only a weak dependence on Mh for a given re- 


If, indeed, our ansatz that matter redistribution 
happens only within is correct, then the total halo 
mass Mh has only a modest influence on the final 
core sizes in dwarf galaxies. Many observations of 
dwarf galaxies in the local universe c an only measure 
dynamical masses ou t to a few kpc (e.g .lDalcanton et all 
20091: lOh et al.ll2(nH iKuzio de Narav et al.ll2006L l2008t 
McConnachiell2012D and so cannot accurately determine 
the halo virial mass. This leads us to Figure [5] where 
we have plotted the range of core sizes as a function 
of stellar mass required to produce a core of that size, 
without adopting a specific coupling efficiency in dark 
matter core formation. Each curve represents the track 
in Vc-Mg space for each halo mass, and is determined by 
the relations between Mh, rh, and Ts (or concentration 
parameter). The endpoints of each curve represent the 
upper and lower limits that we imposed when computing 
the energy requirements from Equation (l5l) . 

In this figure, the efficiency e is equivalent to sliding a 
given dwarf galaxy along the Mg axis. To illustrate this, 
we show the local gr oup dwarf Fornax, w hose stellar 
mass was taken from iMcConnachiel (I2012D and whose 
dark m atter core size was measured bv lAmorisco et ^ 
(1201,31 ). Fornax lies significantly below the predicted 
core size for its stellar mass at 100% efficiency, but is 
consistent with our core predictions if only a few per cent 
of the total SNe budget contributed to the formation 
of the dark matter core in the host halo. Fornax’s host 
halo could range anywhere from 10®-10^^ Mq while e 
could range anywhere from 0.1-1%. 


5. DISCUSSION AND CONCLUSIONS 

In this paper we presented a new estimate for the en¬ 
ergy requirement for galaxies to form dark matter cores 
in their host halos. We assume former cusp material re¬ 
mains within a characteristic radius, r^, at which both 
the density and mass profiles of the cored dark matter 
halo and the original cusped halo match. This signifi¬ 
cantly reduces the amount of energy requ ired compared 
to the estimate of lPefiarrubia et al.l (|2012l) without hav¬ 
ing to requir e, for example, that co res form only at high 
redshift (e.g. lAmorisco et 311120141 ). Our results suggest 
that dark matter cores are present in most, if not all, 
dark matter halos that ex perienced clustered st ar forma¬ 
tion within their centres ([Maxwell et al.l[2012D . In the 
case of the Local Group dwarf galaxy Fornax, if it lives 
in a 10^° Mq halo then less than 1% of the energy re¬ 
leased by all supernovae over its star formation history 
would be required to build a 1 kpc core. 

Our new calc ulations alleviate another potential prob¬ 
lem raised by iPenarrubia et al.l : the accretion of ad¬ 
ditional high phase-space densi ty collisi o nless particles 
(cuspy satellites) was shown bv iDehn^ (I2005D to tend 
to preserve the steepest density cusp in the merger. How¬ 
ever, if all large halos merging in the collisionless mass 
assembly process are already cored, the dark matter den¬ 
sity profile of the merger remnant should be cored as well. 
While the smallest accretors, common at early times, 
may be cuspy and rebuild a central cusp, subsequent star 
formation feedback should be enough to destroy the dark 
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Maxwell, Wadsley, & Couchman 



Fig. 3.— The estimate of the amount of energy /\E required to convert a dark matter cusp into a core. The solid lines with symbols 
show the energy required to form a core of a given fixed size using the pseudo-isothermal density profile (see text). The solid black line 
shows how A.E scales w ith Mh if rm is allowed to scale as a fixed fraction of rh- The grey area corres ponds to the energy e stimate of 
IPeharrubia et ^ l|2Q12 lb The do tted, dashed, and dot-das hed lines corres pond to the Ms-Mh relations of fBehroozi et ahl (|20iCll b [Kravtso^ 
H2010lb and IMoster et ^ l|201QI) respectively, assuming a IKroup^ 1)20011) IMF. The right axis shows the stellar mass corresponding to 
AF? assuming 100% efficiency. 
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Fig. 4.— The predicted core size as a function of halo mass 
using the three relations discussed in this paper. It is 

extremely difficult to compare these predictions to observed dark 
matter cores in local group dwarf galaxies given the uncertainties 
in the relations, and the difficulty in determining for individual 
galaxies. 

matter cusp in a short tim43- 
With the apparent tension between ACDM and the 
presence of cores in dwarf galaxies relieved, we can now 

^ The work of IMashchenko et al.l ()20Q8I) . for example, showed 
that a core can form in 100 Myr. 


Fig. 5.— The predicted dark matter core as a function of stellar 
mass. The grey lines show the energy curve in Vc-Ms space for a 
given Mh- The top axis shows our assumed conversion from Ms to 
SNe energy. For clarity, we have highlighted the curves for halos 
of 10® M©, 10® M©, 10^® M©, and 10^^ M©. The star shows the 
stellar mass (2 x 10^ M©) of Fornax fr omlMcConnachid 112012) with 
the dark matter core size estimated bv lAmonsc^^^aTTll2^^ . The 
efficiency e to which the SNe feedback couples to the creation of 
a dark matter core can be estimated by sliding the Fornax data 
point to the left. 

focus on making firm predictions of the sizes of dark mat- 
ter cores as a function of stellar mass (Figure [5]). For 
example, simulations of dark matter core formation in 


































The Energetics of Cusp Destruction 


7 


dwarf galaxies can test the efficiency and energy output 
of star formation prescriptions. The resulting density 
profiles and other properties can b e directly compared to 
observations of local dwarf s (e.g. iKuzio de Narav et al.l 
120081: iDalcanton et al.ll2009ll . It is also important to note 
that we have neglected all other forms of star forma¬ 
tion fee dback aside from supernovae fe.g. iHookins et al. 
2011 : ?: lShen et al.ll2010HAgertz et al.ll2013HKeller et al. 
201411 . However, these extra sources of energy can simply 
be approximated by tuning e, the efficiency with which 
the energy is transferred to the dark matter, using the 
energy from supernovae as a baseline. 

The ability for a given halo to form a core will de¬ 
pend on the coupling between the depth of the central 


halo potential and the star formation rate which lim¬ 
its the amount of available feedback energy. The bursty 
star formation at the centres of dwarf galaxies should 
be able to drive e h igh enough to fo rm a sizeable core 
(lMax^relI_et_a]J[2Ql^ at high redshift (iMas ichenko et al.l 
120081 lAmorisco et al.l 12014: iMadau et al.l 1201^ How- 
ever, low mass halos may not be able to form cores if 
their galactic potential is shallow enough that each burst 
of star form ation results in completely efficient gas ex¬ 
pulsion fe.g.lNava rro et allllOOPall Read fc Gilmorell^005l : 
iShen et al ll20l4 IMadau et al.ll201'lh 

The authors wish to thank NSERC for support and 
R. Kuzio De Naray for providing the observed rotation 
profile data. 
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APPENDIX 

A. VELOCITY PROFILE FITS 


Figure El repeats the an alysis shown in Figure [T] for all nine Low Surface Brightness galaxies observed by 
IKuzio de Narav et al.l (|2008ll . 
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Fig. 6.— The velocity and density profiles for all of the lKuzio de Narav et ^ (120081) LSBs. Line styles and symbols are the same as in 
Figure fn 
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Fig. 7. — The uncertainty contours for all of the lKuzio de Narav et Il20n8l ~l LSBs, with black and gray indicating the pseudo-isothermal 
and fixed slope profiles, respectively. 


B. ERROR ESTIMATION FOR BEST FIT VALUES 
We determine the best fit by finding the minimum in the square of the velocity residual, 7Z^: 


N 




2=1 


(Bl) 


where Rj is the ve l ocity residual at each measured radius and is the measured velocity error from 
iKuzio de Narav et alJ (|2008ll . If the density models used in the text were linear in the parameters Ts, a, and S, 
the uncertainties could be derived by using the probability distribution to determine confidence levels. Since this 
is not the case, the probability distribution of TZ^ cannot simply be approximated to that of ■ Instead, we compute 
an uncertainty contour given by: 


min (7^^) 
N-M ' 


(B2) 


where M is the number of free parameters — two for Equations @ and (ISl) and one for Equation (IT|) — since the 
central density is fixed. In Figure |B1 these contours a re shown by the solid black a nd gray lines for Equations (I2|) and 
©, respectively, for the sample of nine galaxies from iKuzio de Narav et all (1200811 . Velocity data points that showed 
large variations from the expected velocity profile were ignored when computing the error contours; these are shown as 
open circles in Figure O The location of the minimum Ti? for each are given by the filled triangles. For most of these 
galaxies, minimizing the residuals gives relatively tight constraints on the scale radius and power index for the density 
profiles. However, UGC 4325 and DDO 64 do not have closed contours, which is unsurprising given that constraining 
Ts requires the velocity profile to begin to turn over — traits that are not exhibited in the velocity data for these 
galaxies. Instead, we can only define acceptable lower limits for the density profile parameters. 






























